Staggered-vorticity correlations in a lightly doped t-J model: a variational approach. 
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We report staggered vorticity correlations of current in the d-wave variational wave function for the 
lightly-doped t-J model. Such correlations are explained from the SU(2) symmetry relating d-wave 
and staggered-flux mean-field phases. The correlation functions computed by the variational Monte 
Carlo method suggest that pairs are formed of holes circulating in opposite directions. 



Flux phases have been proposed as mean-field solu- 
tions to two-dimensional antiferromagnets Q. For the 
undoped Heisenberg antiferromagnet, the staggered-flux 
variational wave function gives a relatively good energy 
, even without the Neel long-range order. However, 
for doped systems staggered-flux phases would normally 
break the translational and the time-reversal symmetries, 
except at the flux value of w. Besides, the physical mean- 
ing of the flux is not transparent. It has been suggested 
that flux is related to the spin-chirality correlations ||, 
but such correlations are very complicated for both ex- 
perimental and theoretical study. 

Remarkably, in the case of doped t-J or Hubbard 
models, there is one more indication of a staggered- 
flux phase. Namely, the current-current correlations 
may show the staggered-flux pattern inherited from 
the mean-field phase. We find such a pattern in the 
Gutzwiller-projected d-wave variational wave function for 
the t-J model. Those correlations may be explained 
by the SU (2) equivalence between the d-wave-pairing 
and staggered-flux phases. Finally, we interpret the 
staggered-flux structure of current correlations as a pair- 
ing between holes of opposite "staggered-vorticity" . 
Consider first the d-wave variational wave function 

k 

where Pq is the projection onto the states with a fixed 
number of electrons and without doubly-occupied sites. 
The coherence factors Uk and Vk are of the BCS form for 
d-wave pairing: 

_ Vk_ _ A fc 
Gk ~ u k ~ 6 + Vti + A fc ' 
£fc = — 2(cosfc a; + cos k y ) — /i, (2) 
Afc = A(cosfca; — cos k y ). 

For simplicity, we consider the case fi = 0. Tuning /i 
improves the energy by only 0.2%, and the current cor- 
relations reported in this paper remain practically the 
same. 

Using the variational Monte Carlo (VMC) method, we 
computed the current-current correlations for the wave 
function (pi) in a finite system (2 holes in the 10 x 10 



lattice with periodic-antiperiodic boundary conditions, 
A = 0.55). The current on a link (ij) is defined as 

jij = —jji = i{ c ai c aj ~ C ajCai)- (3) 

The correlation function exhibits a staggered-flux struc- 
ture (Fig. |l|). Since the current vanishes as the hole con- 
centration x —* 0, we have divided the correlation func- 
tion by x. This behavior of the current correlations is 
observed in the whole range of the gap values A. The 
correlations are smaller at small A, which will be in- 
terpreted as the limit of staggered flux going to zero. 
We introduce vorticity V for any plaquet as a sum of 
the currents around it in the counterclockwise direction: 
V = (ji2 + J23 + 334 + ja)n- The vorticity correlations 
(Fig. |lp) obtained from the current-current correlations 
in Fig. [j]a have alternating sign with a phase shift of w, 
so that the sign of (V(0)V(R)) is 

This staggered-vorticity correlation is a surprising con- 
sequence of the projection, because it is absent in the un- 
projected d-wave wavefunction. In order to understand 
this, we show below that the same wavefunction can be 
written as a projection of a staggered-flux state. This 
has the advantage that properties that are obscure in 
one representation may become obvious in another. For 
example, the staggered-flux wavefunction is an insulator 
with gap nodes. The appearance of superconductivity 
after projection is a surprise, but vorticity and attrac- 
tion between holes appear naturally, as we shall discuss 
below. 

The basic starting point is the SU(2) symmetry in the 
fermion representation of the t-J model. This is well 
understood in the undoped case §, where SU(2) dou- 
blets V>Ti = (frufli) and ipii = (fiu-tfi) represent the 
destruction of spin up and spin down in the subspace 
of one fermion per site. Wen and Lee @| extended this 
symmetry away from half filling by introducing a doublet 
of bosons bi = (&H,&2i)- The physical electron is repre- 
sented as an SU(2) singlet formed out of the fermion and 
boson doublets c a i — -^b^ai- 

At the level of variational wave functions, the con- 
straint of no double occupation is enforced by projecting 
the fermion-boson wave function onto the SU (2)-singlet 
subspace of the Hilbert space. On each site, there are 
only three physical states: spin-up, spin-down and a hole, 
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IT) = / T T |o), |j) = /f|o), |*> = ^(b\ + 4/1/1)10). 
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The projector may be written as 

^( 2 )=II(it)(ti + u)ai 



(4) 



(5) 



It should be applied to a mean-field state with the total 
number of bosons defining the number of doped holes. 
The mean-field Hamiltonian has the form 

H = E [ J ^ U ^j + t(b\u l3 b 3 + h.c.) 



(6) 



where Uij are 5*?/ (2) matrices representing generalized 
hopping amplitudes (mean-field parameters) for nearest- 
neighbor sites i and j. 

In the underdoped region, the mean field solution is 
the staggered-fhrx phase characterized by 



0. 



(7) 



with ciij — (— l) a W + "0')^j/4 forming a staggered-flux pat- 
tern around plaquets of the lattice, as shown in Fig. |^a 
(the overall normalization of £7y is of no importance for 
the wave function). The gauge- invariant variational pa- 
rameter of the staggered-flux ansatz is the flux per pla- 
quet ip — ^2 n a,ij. Even though the ground-state wave 
function breaks time-reversal and translational symme- 
tries, these symmetries are restored after the projection 

Even with hole doping, the fermion bands are exactly 
half-filled, and therefore the number of "no-fermion" sites 
is equal to the number of "two-fermion" sites. This is 
shown in Fig.|^b. Note that both these sites are spin 
singlets and have the right spin quantum number for a 
physical hole. In SU(2) theory, a b\ boson is attached 
to the "no-fermion" site and a &2 boson to the "two- 
fermion" site, and both become physical holes, according 
to Eq. (|J). In the mean field theory, the bosons con- 
dense to the bottom of their respective bands, which 
are located at Qi = and Q2 = (n, it). The pre- 
scription of constructing an SU (2) projected wavefunc- 
tion is as follows. The physical wavefunction is speci- 
fied by the location of the up spins and the holes (cir- 
cled sites in Fig.||b). A given set of holes is partitioned 
into all possible "no-fermion" and "two-fermion" sites, 
denoted by {ri},{r2}. Each partition specifies a con- 
figuration of the staggered flux wavefunction given by 
a product of two Slater determinants (for spin-up and 
spin-down spinons). To this we multiply the phase fac- 
tor exp[iJ2{ ri } {ra}^ 1 ' r i + Q2 • r 2 )] to represent Bose 



condensation of 61, 62 and sum over all partitions. An ad- 
ditional sign depending on the ordering of sites is needed 
to preserve the antisymmetry of the wavefunction. This 
prescription realizes the projection Psu(2) of a mean-field 
staggered-flux state to the SU (2) singlet subspace. 
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(b) (c) 
FIG. 1. (a) Current-current correlations for 2 holes in the 
10 x 10 lattice at A = 0.55 and fi = 0. Boundary conditions 
are periodic in one and antiperiodic in the other direction (the 
data are averaged over the two orientations). The number on 
a link is the correlation of the current on this link and of the 
current on the circled link divided by hole density. The ar- 
rows point in the direction of the positive correlations of the 
current, (b) The same data in the form of vorticity correla- 
tions. The number on a plaquet is the vorticity correlation 
divided by x with the crossed plaquet. (c) Same as (b) for 10 
holes in 10 x 10 lattice. 

Next we prove the equivalence of the SU{2) projected 
staggered flux phase with the pairing state (|l|). Such 
an equivalence is known in the undoped case we 
extend it to doped systems. We note that the mean- 
field ground state of the staggered flux phase has the 
form |$„) = |#/> ® |$ 6 >, where |$ 6 ) = e^ B l+^ B l\0) 

with {Bi,B 2 ) = (Ei^EiC-l)^), and l $ /> is thc 
fermion state in the staggered flux phase Eqs. @,(0). 
Complex numbers b\ and 62 parametrize the Bose con- 
densate where b\ y 2 condense into their band bottoms: 
B a \$b) = b a \®b}- The SU(2) projection selects the sub- 



2 



space with equal numbers of Bi and B 2 bosons, and 
therefore states corresponding to different choices of b% 
and &2 differ after projection only by a one-parameter 
transformation exp(AiV/j), where A is a complex param- 
eter, Nh is the number of holes. If we fix the average 
hole concentration, the wave function is defined unam- 
biguously, up to an unimportant gauge exp(AiV^) with 
ReA = 0. 
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(a) (b) 
FIG. 2. (a) Staggered-flux phase. Links with arrows cor- 
respond to ciij = if/4 in the direction of the arrow, (b) A 
typical configuration of the half-filled fermion state. Arrows 
denote fermions. Circled sites are physical holes which are 
spin singlets made up of either empty or two-fermion sites. 61 
and &2 bosons are assigned to these respective sites. 



The staggered-flux mean-field state is related to the 
d-wave state by a SU(2) rotation Wi = exp[(— l) J; ?ri]: 



w iUijW ],[ ± l /2 ±A / 2 ) 



(8) 



where the sign of A is opposite for vertical and horisontal 
links, and A related to tp by 
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2"' 



(9) 



After the SU(2) rotation, the mean field state has the 
same form, except |<E>/) is replaced by |4>y), the fermion d- 
wave state of the BCS form (Q), and the bosonic parame- 
ters b a are rotated: b' x = (&i+i&2)/v2, b' 2 = (ib~i+b 2 )/\/2. 
Since the two mean-field states are related by a SU(2) 
gauge transformation, they lead to the same physi- 
cal state after the SU(2) projection: Psu(2)\^mean) — 
Psu{2)Wmean)- The freedom in the choice of the bosonic 
parameters b\ and b 2 established in the staggered-flux 
gauge allows us to set b' 2 = 0, then the SU(2) projec- 
tion becomes equivalent to the conventional Gutzwillcr 
projection. This proves that the SU (2) wave function is 
identical to the conventionally-projected wave function. 
Of course, the above proof equally applies to the systems 
with fixed numbers of holes which we use in our VMC 
calculation. In this case, the SU (2) projected wave func- 
tion is identical to that given by Eq. (Q). 

Note that in the proof of the equivalence of the two 
wave functions we use a^ 1 = (and, therefore, fx = ). 
A finite fi in the pairing gauge corresponds to a nonzero 
value of the Lagrange multiplier <zf (chemical potential). 



In the staggered-flux gauge, this translates into an on-site 
pairing term 



(10) 



with some site-dependent phases 0j. This term only 
slightly affects the properties of the projected wave func- 
tion, and we neglect it in the further discussion. 

In the ground state of the staggered-flux mean-field 
Hamiltonian [eq.@ with Uij given by (0)], the b\ and 
62 bosons attract each other. This can be understood 
from the correlation function for the "excess density of 
fermions" : 

((1 - n f (i)) (I - n f (j))) = ((l - flj m ) (l - fy M )) 



l(4/oi)l 2 <o 



(If) 



at the mean-field level. This means that around a "no- 
fermion" hole (61 hole) with 1— 71/ (i) > there is a region 
of an increased probability to find a "two-fermion" hole 
(62 hole) with I — 71/ (j) < 0, and vice versa. The mean- 
field Green's function G(i,j) = {f ai faj) decays as Br 2 
at large distances R = \i — j\. This is a consequence of 
the nodes k = (±7r/2, ±7r/2) in the mean-field spectrum 
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(12) 



Thus at the mean-field level, the attraction of the two 
species of holes leads to a R~ A decay of the density- 
density correlations. After the projection, the attraction 
becomes much weaker, but still survives, as we shall see 
from our VMC computations. This attraction between 
holes was observed earlier || , but was difficult to explain 
in the d-wave gauge. 

In Fig. H we plot the correlations of the hole density 
Uh(i) = 1 — CaiCai and of the vorticity V at A = 0.55 
for two holes in a 18 x 18 lattice as functions of dis- 
tance. We observe power-law decay of both correlations, 
(n h (0)n h (R)) cx (V(O)V(R)) cx with the equal ex- 

ponents a ~ 1.2. 

The reduction of the exponent a from its mean-field 
value a = 4 is due to the projection. For the case 
of two holes ("zero-doping limit"), this reduction is so 
strong that a < 2 and, as a consequence of the sum rule 
^2ji{ n h(0)nh(R)) — nh , the two holes become unbound 
as the size of the system increases. 

The proportionalty of the vorticity and density corre- 
lations, together with the sign of the vorticity correla- 
tions, suggests that the two holes may be thought of as 
carrying opposite staggered vorticity. This can be jus- 
tified if we fix the gauge and describe the wavefunction 
as the staggered-flux phase. At the mean-field level, this 
breaks the time-reversal symmetry, and the fermions in 
the filled single-particle states carry a non-zero staggered- 
vorticity with (jl) ^ 0, where jl = i(flj a j - ftjai) is 
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the (unphysical) fermion current. The SU(2) projection 
restores the time-reversal symmetry, and for the pro- 
jected wave function (jij) = 0, for the physical current 
jij. In the staggered-flux gauge (0), at the mean- field 
level, jij ~ j{j for b\ holes and jij ~ I0r ^2 holes. 
Thus (jij) = for the physical current is a result of 
the balance between the opposite staggered-vorticity of 
61 and b-2 holes. The attraction between the two SU(2) 
species of holes then implies attraction between holes cir- 
culating in the opposite directions. For a finite system 
with two holes, the holes are always of opposite SU(2) 
types, which resultins in the proportionality between the 
vorticity and density correlations. 

The above picture can be summarized in the expression 
V cx pi—p2, where V is the staggered-vorticity (— ) R V(R) 
for the physical current jij, and p\$ are the densities of 
the two bosons 61,2- The correlations of the staggered- 
vorticity and of pi — p 2 are related (V(R)V(0)) cx 
{ Pl (R)pi(0)) + {Pi (R)P2 (0)}-2{px(R)p 2 (0)). When there 
are only two holes (pi(R)pi(0)) — (p%{R) P2(0)) = 0, we 
find that 

(V(R)V(0)) cx -(p(R)p(0)) (13) 

where p = pi + P2 is the total density of the holes. The 
minus sign explains the it phase shift seen in Fig. 1. 




t 10 100 

FT2 



FIG. 3. Hole density and staggered-vorticity correlations 
for 2 holes and staggered-vorticity correlations for 16 holes 
in the 18 x 18 lattice. Boundary conditions are periodic in 
one and antiperiodic in the other direction, A = 0.55. The 
correlation functions are divided by the density of holes x 
and plotted as a function of the squared distance. The data 
are obtained as a result of averaging over 2 ■ 10 4 samples for 
2/18 x 18 system and 2 ■ 10 3 samples for 16/18 x 18 system. 

If we assume pairing between holes of opposite stag- 
gered vorticity, we may interpret the vorticity correla- 
tions as the hole correlations within one pair which allow 
us to determine the strength of the pairing correlation 
(even for finite density of holes) and measure the size of 



pairs. At finite density of holes, the correlation of the 
staggered-vorticity decay faster. In Fig. 4b we present 
the correlation for 16 holes in 18 x 18 lattice (nearly 5% 
doping) at A = 0.55. We find a as 2.2. The increase of a 
above 2 may be interpreted as a formation of bound hole 
pairs (and hence the onset of superconductivity). The 
small value of a — 2 after the projection implies that the 
pairs are bound very loosely and their size is large. In 
this case, the nearest-neighbor pairing amplitude (c^Cji) 
is small, which may account for numerical conclusions 
about the absence of superconductivity at J/t < 0.5 

While the staggered vorticity correlation is found for 
a superconducting wavefunction, we speculate that the 
phase coherence of the bosons may not be crucial and 
that such correlation may survive above T c in the pseu- 
dogap state, and indeed serves as a signature of that 
state. Unfortunately, detection of this correlation may 
be difficult. We estimate that the fluctuating current 
generates a fluctuating staggered magnetic field of order 
40 G. This field will contribute to the relaxation of the 
Y nuclei, which are ideally sited above the center of the 
plaquets. However, we do not have dynamical informa- 
tion at present and it is difficult to predict the magnitude 
of this orbital relaxation. Another possibility is to freeze 
in some staggered field pattern around an impurity site. 
If some Y is replaced by an impurity with spin S which 
strongly couples to the Cu-0 plane, we may use a mag- 
netic field to align the spin and therefore create a static 
orbital current via the L • S term. This local orbital cur- 
rent may then generate a static staggered pattern around 
it, which may be detected by shifts in the Y NMR line. 
Perhaps a more promising way to test our prediction is 
to look for this effect in exact diagonalization of small 
systems. 
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